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Using Effective Field Theory (EFT) methods we present a Lagrangian formalism which describes 
the dynamics of non-relativistic extended objects coupled to gravity. The formalism is relevant to 
understanding the gravitational radiation power spectra emitted by binary star systems, an impor- 
tant class of candidate signals for gravitational wave observatories such as LIGO or VIRGO. The 
EFT allows for a clean separation of the three relevant scales: rs, the size of the compact objects, 
r the orbital radius and r/v, the wavelength of the physical radiation (where the velocity v is the 
expansion parameter). In the EFT radiation is systematically included in the v expansion without 
need to separate integrals into near zones and radiation zones. Using the EFT, we show that the 
renormalization of ultraviolet divergences which arise at in post-Newtonian (FN) calculations 
requires the presence of two non-minimal worldline gravitational couplings linear in the Ricci cur- 
vature. However, these operators can be removed by a redefinition of the metric tensor, so that 
the divergences at arising at have no physically observable effect. Because in the EFT finite 
size features are encoded in the coefficients of non- minimal couplings, this implies a simple proof of 
the decoupling of internal structure for spinless objects to at least order . Neglecting absorptive 
effects, we find that the power counting rules of the EFT indicate that the next set of short distance 
operators, which are quadratic in the curvature and are associated with tidal deformations, do not 
play a role until order v^". These operators, which encapsulate finite size properties of the sources, 
have coefficients that can be fixed by a matching calculation. By including the most general set of 
such operators, the EFT allows one to work within a point particle theory to arbitrary orders in v. 



I. INTRODUCTION 



The experimental progam in gravitational wave detec- 
tion that is currently under way 0, Q has brought re- 
newed attention to the problem of obtaining high accu- 
racy predictions for the evolution and gravitational radi- 
ation power spectra of binary star systems with neutron 
star or black hole constituents. The conventional ap- 
proach to calculating the initial inspiral of the system as 
it slowly loses energy to gravitational radiation employs 
the post-Newtonian (PN) approximation to general rel- 
ativity. Essentially, this formalism consists of systemati- 
cally solving the Einstein equations with non-relativistic 
(NR) sources as a power series expansion in w < 1, where 
w is a typical three-velocity of the system under consid- 
eration 3, 4J. For binary systems that can be detected 
by LIGO, which have binary constituents with masses in 
the range ^ (1 — 10) tuq, the inspiral phase spans roughly 
0.1 < w < 0.4, corresponding to orbital separations from 
about 100 to a few times {m/mo) km. 

A detector such as LIGO is particularly sensitive to 
the time varying phase A(/)(t) of the gravitational waves 
emitted during the evolution of the binary system . In 
the PN formalism, this observable can be obtained by 
calculating the mechanical energy E as well as the total 
power P emitted in gravitational radiation as functions 
of the orbital parameters of the binary system. Taking 
for illustration the extreme mass ratio limit and using 
conservation of energy, dE/dt = — P, for circular orbits 
the gravitational wave phase seen by the detector is then 



given by 

Ac^it) = fdMr) ^ -1- r dv'v'^^§l^, (1) 
Jto GnM P[v') 

where u}{t) is the wave's frequency and M the mass of 
the heavy component star. Because a typical inspiral 
event will sweep a large number of radians of phase as it 
scans the detector frequency band (e.g, for LIGO, with 
a frequency range 1 — 10^ Hz, the number of cycles for a 
neutron star/neutron star binary is of order 10^), to com- 
pare the predictions of GR with the experimental data 
requires computing the observable A(/)(i) to extremely 
high order in the velocity expansion. For instance, for bi- 
nary neutron star inspirals in the LIGO band, one must 
know A0(i) to 0{v^) (or "3PN") beyond the leading or- 
der quadrupole radiation results, and higher order terms 
in V may be necessary for tracking the phase variation 
of the gravitational waves emitted by more massive ob- 
jects |3- 

Several independent research groups have applied vari- 
ants of the PN formalism to the computation of high or- 
der corrections to the gravitational wave observablcs dur- 
ing the binary inspiral phase 6, 7. ^ lol Il0lllli |. Heroic 
calculational efforts have yielded in the last few years the 
analytic results necessary to fully interpret the LIGO 
data [lOl llr|. While the independent calculations are in 
agreement including terms up to (or "2.5PN") beyond 
the predictions of Newtonian gravity, the corrections 
do not give well defined results. At this order, there 
seems to be an ambiguity in the PN expansions for the 
binary system, stemming from scheme dependence in the 
regularization of singular integrals that arise in the per- 
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turbative series. Recently, refs. [T^ IT^ have proposed a 
resolution of this ambiguity purely within the framework 
of the PN expansion. 

The singular integrals that cause trouble at high orders 
in V can be seen to arise from treating the binary star 
constituents as point particle sources coupled to general 
relativity. It is then clear that these singularities are sim- 
ply the usual ultraviolet (UV) divergences that appear in 
any field theory coupled to point sources. Their presence 
is merely a sign that the formalism has limited predictive 
power and needs to be supplemented by a more complete 
model that accounts for the fact that the sources are not 
really pointlike. While in principle one could obtain fi- 
nite, unambiguous results for experimental observables 
in this more complete theory, the price would be much 
greater computational complexity. Even then, it would 
be difficult to disentangle the effects that depend upon 
short distance physics from those that are solely a con- 
sequence of the gravitational dynamics. Since the goal 
of the gravitational wave experiments is both to probe 
the internal structure of strongly gravitating astrophys- 
ical systems as well as to test the dynamics of general 
relativity (or possible deviations from it), it is critical to 
have a computational approach that can cleanly separate 
these two aspects of the binary problem. 

In this paper, we propose a Lagrangian framework for 
systematically calculating to any order in the PN expan- 
sion. Within our theory the "ambiguities" that plague 
the conventional PN techniques at order and beyond 
have a natural place in the renormalization program. Our 
formalism is based on effective field theory (EFT) rea- 
soning, and it is motivated by the EFT approach to NR 
bound state problems in quantum field theories such as 
QED and QCD jTHH. Rather than trying to resolve 
the point particle singularities by resorting to a specific 
model of the short distance physics, in an EFT frame- 
work 

m we systematically parametrize our ignorance of 
this structure by including in the effective point particle 
Lagrangian the most general set of operators consistent 
with the symmetries of the well understood, long wave- 
length physics. In our case, this long wavelength physics 
is general relativity, and the symmetry that constrains 
the dynamics is just general coordinate invariance. 

The operators in the point particle action have coeffi- 
cients which encapsulate the properties of the structure 
of the extended sources. UV divergences encountered in 
calculations in the EFT can then be dealt with using 
any convenient regulator, and renormalized by adjust- 
ing the operator coefRcients^ . Given a model of stellar 
structure, the precise value of these coefficients can be 



obtained by a UV matching calculation that we describe 
below in Sec. IIVI Note that because these operators are 
intrinsic to the structure of the sources, the matching pro- 
cedure can be carried out for a single source in isolation, 
independently of the complicated binary star dynamics. 
In this way the EFT manages to disentangle the model 
dependent aspects of the gravitational wave signals from 
the features that arise as unambiguous predictions of gen- 
eral relativity. 

An EFT also has the advantage of having manifest 
power counting in the expansion parameter of the the- 
ory, in our case the velocity v. This means that in the 
EFT it is possible to calculate, using simple scaling argu- 
ments, the order in v which a given term in the perturba- 
tive series (i.e, a typical Feynman diagram obtained from 
the vertices of the effective Lagrangian or an operator in 
the effective theory) first contributes to a given physical 
observable. In particular, one can determine ahead of 
any calculation the number of terms in the point particle 
Lagrangian that will appear at any order in v. 

As an example, using our formalism we will sec that 
for binaries with non-spinning black hole or neutron star 
constituents, the proper time term is sufficient to de- 
scribe the physics up to v^. The apearance of logarith- 
mic divergences in the PN expansion that first appear 
at order indicates that two new non-minimal terms 
whose form is dictated by general covariance must be in- 
cluded at this order to consistently renormalize the point 
particle theory. However, the coefficients of these oper- 
ators can be shifted arbitrarily by making field redefini- 
tions of the metric, which indicates that the operators 
cannot contribute to observables. This gives a clean ar- 
gument for the absence of finite size effects in the 
predictions for gravitational wave templates. Ignoring 
the effects of absorption^, we also find that the next set 
of point particle operators, including for instance oper- 
ators quadratic in the Riemann curvature, do not con- 
tribute until order v^'^. These operators, which cannot 
be removed by field redefinitions, are associated with the 
leading effects of tidal forces exerted on an isolated neu- 
tron star or black hole by an external gravitational field. 
The fact that such tidal operators do not arise until v^^ 
is consistent with purely Newtonian estimates which also 
conclude that tidal forces are w^" effects'^. 

Finally, an EFT has the advantage that it is possi- 
ble to use renormalization group (RG) methods to un- 
derstand in a simple way logarithms of scale ratios that 
may arise in the perturbative series for a given observ- 
able. For example, in the NR expansion of the binary 
star gravitational wave observables, one encounters log- 



In this paper we employ dimensional regularization, and modified 
minimal subtraction (the MS scheme). This has the advantage 
of not introducing spurious length scales that c ould spoil dimen- 
sional analysis arguments within the EFT. See |13| for a discus- 
sion of dimensional regularization in the context of conventional 
PN formalisms 



^ In the case of extreme mass ratio binaries where the large mass 
component is a black hole, it is kno wn that absorption into the 
black hole first appears at order llTl . We will show how ab- 
sorptive effects are incorporated into the EFT in a future pa per . 

^ We thank Kip Thorne for pointing this out to us. See also |1S| 
for a recent review. 
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arithmic terms of the form t;"(lna)™, where a = rg/r is 
the ratio of a distance scale that characterizes the size 
of the binary constituents and r is a distance of order 
the orbital radius (note that for black hole binaries, the 
virial theorem relates a ~ v^). We will show that in our 
EFT, the coefficients of operators in the effective point 
particle action undergo non-trivial RG flows, so that the 
logarithms of a (tidal logarithms) can be obtained by RG 
running from a renormalization scale fi ^ 1 jr^ where the 
point particle theory is initially defined down to a scale 
[I ^ Xjr where the binary dynamics takes place. More- 
over, the complete structure of the RG flows is encoded 
in the black hole metric. 

As an explicit example of the RG method, we will give 
a simplified derivation of the terms at O(w^) in the two- 
body interaction Lagrangian (or equivalently the binary 
gravitational binding energy) that are enhanced by In a. 
We then show explicitly that the effect of these loga- 
rithms do not contribute to physical observables. This 
result must follow since, as we already mentioned above, 
the worldline operators can be removed by redefining 
the metric. However, this calculation is still useful as 
a template for determining the logarithmically enhanced 
terms at higher orders in v that stem from true finite size 
effects. 

The outline of the paper is as follows. In the next 
section we formulate our EFT. The starting point is a 
theory of relativistic point particles coupled to gravity. 
This theory is appropriate for discussing the dynamics of 
extended sources at inverse momentum transfers larger 
than their size r<,. We then describe the general proce- 
dure for matching this theory at the mass scale* /i ^/r, 
where the binary orbital dynamics takes place, onto an 
EFT that is optimized for the NR limit (NRGR). The re- 
sulting EFT has manifest power counting in the typical 
velocity v of the sources, and governs the dynamics of the 
long wavelength "radiation" gravitons that propagate out 
towards the gravitational wave detector. To illustrate our 
formalism, in this section we also give two classic exam- 
ples of how the matching procedure and the EFT velocity 
power counting rules are applied to reproduce the leading 
non-trivial consequences of general relativity in the con- 
text of the NR two-body problem: the Lagrangian that 
describes the gravitational interactions of two NR objects 
at 0{v'^) beyond Newtonian gravity Q, and the interac- 
tions of gravitational radiation with the binary star at 
leading order in the velocity. 

In Sec, mil we return to the relativistic point particle 
effective Lagrangian, and show how all point particle di- 
vergences can be absorbed into worldline operators. In 
Sec. IIVI we discuss how, given a more complete stellar 
model one can determine the coefficients of operators in 
the point particle EFT at the scale ^ ~ l/r^. Finally, in 
Sec. we show how the RG running of the EFT coeffi- 



We use units Ti = c = 1. 



cients at order can be exploited to yield a simplified 
calculation of the terms in the NR expansion of the two- 
body gravitational Lagrangian that scale as In a. This 
calculation serves as a prototype for the RG running of 
worldline operators. We have also included an appendix 
which illustrates how the Feynman rules of the EFT are 
derived and used. 



II. THE EFFECTIVE THEORY 

A. Kinematics 

To construct an EFT that has manifest velocity power 
counting, we must first carefully analyze the kinematic 
configurations that arise in the binary problem. Con- 
sider a binary system composed of slowly moving neu- 
tron stars or black holes. Let m characterize the typical 
mass of the constituents (taking them, in this discussion, 
to be roughly the same size), r the typical separation 
between them and v their relative three-velocity. In a 
system whose evolution is determined to leading order 
by a 1/r Newtonian interaction, these parameters are 
not independent, since the virial theorem relates 



The EFT for the binary system will contain degrees 
of freedom representing the point particle binary con- 
stituents coupled to gravitons. While the typical particle 
momenta are of order {E ~ tow^,p ~ mv), the gravi- 
tons appearing in a generic Feynman diagram have mo- 
menta that can be divided into two classes. Gravitons 
with momentum scaling as (A;° ~ t;/r, k ~ 1/''), me- 
diate the forces responsible for binding the 2-body sys- 
tem. These potential gravitons can never go on shell and 
thus will not appear as propagating degrees of freedom. 
On the other hand, radiation gravitons with momentum 
{k'^ ~ w/r, k ~ v/r) can appear on shell and must be 
kept in the EFT to reproduce the correct long distance 
physics. 

Note that the interaction of an NR particle with a sin- 
gle potential or radiation graviton causes the particle to 
recoil by a fractional amount |k|/|p| ~ h/L <^ 1 , where 
L ^ mvr is the typical orbital angular momentum of the 
system. Thus as far as the graviton dynamics is con- 
cerned, the NR particles can be treated as background 
non-dynamical sources. Consequently, in this regard the 
EFT that describes the binary configuration has more 
in common with heavy particle effective field theories 
(eg., HQET [l^l), in which light modes interact with 
a static, non-propagating source than a theory such as 
NRQCD (1^, J^J, where the quark and gluon degrees of 
freedom have momentum components which differ only 
by powers of the velocity. However, unlike HQET our 
EFT is not an expansion in inverse powers of the mass. 

As we will see, treating the stellar objects as non- 
propagating sources of gravitons also resolves a long 
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standing formal problem in the treatment of general rela- 
tivity as an EFT |^ . For astrophysical systems with 
masses m ^ mp;, one would naively conclude from the 
conventional Feynman rules of the EFT that the expan- 
sion parameter is m/mpi ^ 1 and therefore the energy 
expansion implicit in the interpretation of gravity as an 
EFT seems to be spoiled |23|- In fact this obstruction to 
the EFT interpretation of general relativity is only an ar- 
tifact of treating the point sources as dynamical degrees 
of freedom, and disappears when the theory is formulated 
in terms of worldlines interacting with gravitons. Once 
this is done, the power counting in v becomes manifest 
at the level of the action. 



B. Point particle Effective Theory 

Given these remarks, the starting point of our EFT 
formulation consists of a theory of relativistic point par- 
ticles coupled to gravity 

S = Seh + Spp, (3) 
where the Einstein- Hilbert term 

Seh = -2mli J d^x^R{x), (4) 
describes the graviton dynamics^ and 

^pp = ~ X! / ^■^'^ + X! '^^^ / dTaR{Xa) 

+ J2 ^ / dTaR^,{Xa)x^aK + ' ' ' (5) 

determines the motion of the 2-body system (a = 1,2 
labels the different particles). In this equation, dxa — 

g f^iydxadx^ is the proper time along the worldline 
of the a-th particle. We will ignore in this paper addi- 
tional degrees of freedom that describe the spin or any 
additional multipole moments carried by each particle. 
This is justified for an extended object in which the en- 
ergy cost of exciting internal degrees of freedom is large 
compared to the typical frequency of the gravitons which 
it interacts with. Thus strictly speaking the formalism 
presented here can be used only to describe the dynamics 
of spinless black holes or spherically symmetric compact 
objects. We will save the problem of including multipole 
moments in the EFT for future work. 

The first term in Spp generates geodesic motion about 
the metric gp,^. Besides this term, we have also explicitly 
shown the first two of an infinite set of possible non- 
minimal couplings of the point objects to the spacetime 



^ Our conventions are R^y = 9ar"^ — cJi/F"^ and signature 

(+,-,-,-). ml^ = l/327rGjv 



metric. Since these non-minimal couplings cause devia- 
tions from pure geodesic motion, one would expect them 
to be associated with the effects of finite spatial extent 
of the sources. Thus by including the most general set 
of such operators, we systematically take into account all 
possible corrections due to the finite size of the object. 
As we discuss later, the operators linear in the curvature 
appearing in Eq. Q will not contribute to any physical 
observable^. This should be contrasted with higher di- 
mension operators in the point particle action involving 
higher powers of the curvature tensor, which encode tidal 
deformations of the binary constituents. We will discuss 
in a Sec. II VI how one goes about determining the precise 
relation between the coefficients of these operators and 
the microscopic physics (ie, the stellar structure model) 
that determines the internal structure. A discussion of 
the laws of motion for extended objects in general rela- 
tivity can be found in 25]. 

A gravitational wave detector such as LIGO measures 
indirectly the power emitted in gravitational radiation 
from the binary system (see Eq. |^). Given the ac- 
tion Eq. (Q, it is in principle possible to calculate this 
quantity by expanding the metric around flat space, 
.9/11^ = ^/i!^ + ^At!v, and integrating out the graviton field 
hp^ to obtain an effective action for the particle coordi- 
nates alone 

e-xjp[iSeff[xa]]= J Vhfj^^expliSEH + iSpp]. (6) 

The effective action Seff[xa] has a real part which gen- 
erates the coupled equations of motion for the 2-body 
system, and consequently its mechanical binding energy. 
It also has an imaginary part that measures the total 
number of gravitons emitted by a fixed two-particle con- 
figuration {x^} over an arbitarily large time T ~f oo 

llm^,,,!..]^/.™-^, (7) 

where dT is the differential rate for graviton emission 
from the binary system. From this quantity we obtain 
the classical power spectrum dP = EdT, and therefore 
the gravitational wave phase seen by the detector. 

In principle, one could directly evaluate Eq. ((HJ using 
the Lorentz covariant Feynman rules generated by the 
Einstein-Hilbert Lagrangian (the Feynman diagram ex- 
pansion for perturbative gravity is reviewed, for example, 
in [23, H^), taking the NR limit at the level of each am- 
plitude. However, the perturbative series generated in 
this way is not optimal for taking the limit v 1. For 
example, the one-graviton exchange term in iSeff[xa\, 




These operators have also been discussed in Ref. |24| 



5 



contains an infinite number of terms from the NR ex- 
pansion of the proper times dr^ = r]^,^dx^dx'^ and the 
Feynman propagator Dp{x). While for a term Hke this 
it is easy to do the expansion in smaU velocity, in a more 
complicated diagram with multiple internal graviton lines 
it becomes more cumbersome to keep track of all the nec- 
essary terms at a given order in v. Furthermore suppose 
one is interested, because of the limited sensitivity of the 
detector, in computing an observable only to some fixed 
order in v. Then how many diagrams do we need to keep? 
It is not possible to answer this question with the Feyn- 
man rules generated by Eq. (O as it stands. Likewise, 
given the non-minimal tidal couplings in Eq. ^ it is not 
possible to determine at what order in v they may first 
contribute. Since from the point of view of the point par- 
ticle theory, Eq. the NR dynamics occurs in the far 
infrared, to formulate an EFT that is tailored to the limit 
<C 1, we must integrate out all field degrees of freedom 
with wavelengths shorter than the orbital distance scale 
r. We now turn to a formulation of this EFT. 



C. NRGR 

The EFT that systematically describes the NR two- 
body problem, which we call NRGR, has manifest power 
counting in the velocity v. To construct NRGR, we pro- 
ceed in several steps. We work in a preferred set of frames 
in which the relative velocity is small, allowing us to ex- 
pand in V. Then in order to have interaction vertices with 
homogeneous velocity scaling, we simply expand Spp in 
powers of the particle three-velocities. Ignoring for now 
the effects of the higher-dimension tidal operators, this 
leads to a worldline Lagrangian 



1 2 
— V 

2 " 



1 . 



hoiVai — ^/iooVq 



(9) 



where h^o, hoi, hij are evaluated on the point particle 
worldline {x'^ ,x.a{x^)) (we use the Euclidean metric Sij 
to raise/lower spatial indices i,j = 1,2,3). 

The propagator for the field h^,j appearing here is still 
fully relativistic, so it does not distinguish between po- 
tential gravitons, which have no right to appear as de- 
grees of freedom in NRGR, and radiation gravitons. Be- 
cause of this, the Feynman rules still do not scale ho- 
mogeneously with V. To deal with this problem, it is 
convenient to decompose the graviton as 

(10) 



h^i^^ix) = h^^{x) + H^^{x), 
where i/^j^ represents the potential gravitons, with 



in, 

r 



(11) 



and h^i, describes a long-wavelength radiation field 

- V - 

dah^,u ^ -hp,^. (12) 



Actually, it is better to further decompose H^^ by remov- 
ing from it the large momentum fluctuations k ~ 1/r. To 
do this we simply work with the Fourier transformed field 
la i?kM.(2:°) 



Hp^{x) 



[ e'^'-^Hk^^ix"), 



(13) 



where is shorthand for J c?'^k/(27r)'^. The advantage 
of this redefinition is that now derivatives acting on any 
field in the EFT scale in the same way, 9^ ^ v/r, so 
it is easy to count powers of v coming from derivative 
interactions. 

The NRGR lagrangian can then be derived by com- 
puting the functional integral 

eyip[iSNRGR[xa,h]] = J VHpi, eyip[iS[h+H, Xa]+iSGF], 

(14) 

where Sgf is a suitable gauge fixing term. We have not 
included ghost terms because as we will see, they are not 
needed in any NRGR computation. Eq. (|14|l indicates 
that as far as the potential modes i/^j/ are concerned 
hp,y is just a slowly varying background field. To preserve 
gauge invariance of the NRGR action, we will chose Sgf 
to be invariant under general coordinate transformations 
of the background metric g^v{x) — i]^^ + h^^{x). To be 
definite, we take the gauge 



with r„ 



^GF 



''PI 



(15) 



i_D^_ff^, where Z?^ is the covari- 



ant derivative derived from the background metric 5^^. 
In this gauge, the 0{H^) terms in the action are (af- 
ter performing a rescaling of i/^^ to obtain a canonically 
normalized kinetic term) 



Chi 



1 



k2 

2 



(16) 



where = H^^^. The terms in the second line of this 
equation are suppressed relative to the first line by a 
power of , and are treated perturbatively, as operator 
insertions, in correlation functions. Thus, the propagator 
for iJk^i^ is 



where P, 



(17) 
with 



d the spacetime dimension. If we assign the scaling x^ ~ 
r/u, k ~ 1 /r, we learn from this that a potential graviton 
scales as H]f^py ^ v^l'^r'^ . Likewise, because in mometum 
space the radiation graviton propagator scales as 1/A:^ ~ 
/v^, the coordinate space radiation graviton field is 



taken to scale as h^. 



v/r. 
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(p",p) 



FIG. 1: The incoming potential graviton with momentum p 
absorbs a radiation graviton with momentum k. 



k 






m/mpi 


1/r 


2 1/2 


v/r 





TABLE I: NRGR power counting rules. 



Just as we had to expand in terms of momen- 
tum labels to obtain manifest power counting it is 
also necessary to perform a multipole expansion of the 
radiation field h^^, whenever it couples to particle world- 
lines or to potential modes ^2^. To see the necessity for 
this, consider the absorbtion of a radiation graviton from 
a potential graviton as shown in Fig. ^ The outgoing 
potential graviton has a propagator of the form 



1 



(p + k)2 



(18) 



which does not scale homogenously in v since p ~ 0{l/r) 
and k 0{v/r). The propagator must be expanded in 
powers of |k|/|p|, which is achieved by multipole expand- 
ing the interactions of the radiation graviton field with 
the potentials at the level of the action, 



/i^^(a;°,x) = /i^^(x°,X) -|-(5xi(9jft,^^(a;°,X) 
-f ]^5:x.i5Tijdidjh^^{x^ , X) -I- • 



(19) 



where i5x = x — X, and X is some reference point in the 
vicinity of the point particle ensemble, for instance the 
center of mass. Likewise, consider the amplitude for the 
point particle system to absorb a 00-polarized radiation 
graviton with 4-momentum k 



2 ^ 



mpi J 



— ik-: 



(20) 



A graviton emitted during the binary inspiral phase has 
4-momentum k {v/r, v/r), so if the coordinates x^ are 
measured relative to, for instance, the center of mass, 
then we have x.a/x'^ ^ v and consequently this amplitude 
contains an infinite number of terms each with a different 
velocity scaling. This situation is also remedied by the 
multipole expansion in Eq. (|19|) . 

These field redefinitions are enough to recast the orig- 
inal point particle theory of Eq. (PJ into a form that is 
better suited to a small velocity expansion. The power 
counting rules in Table^allow one to determine the order 
in the velocity at which any diagram in the calculation of 



the functional integral in Eq. H14() , as well as in Feynman 
diagrams computed in NRGR. However the simplest ex- 
ample of the power counting rules, the counting of the 
particle NR kinetic term, reveals a slight subtlety that 
we must deal with. Given that dx'^ ~ r/w, we have 



dx°r 



2 

—mv 

V 



(21) 



Thus this term scales as the orbital angular momentum 
L ^ 1. This is not the only term that scales in this way. 
The exchange of a single Hqq potential graviton between 
two particles, responsible for the Newton force, involves 
the computation of the diagram of Fig.|2K, which is given 
by 



Fig.ISi = J dhdt25{h~h) j^^e 



i / dt- 



GN'mim2 

'|xi(t)-X2(i)|' 

{Gn = l/32nmpi). This scales as 

[{dt){d^k){m/mpi)H]'^ - mvr, 



-ik-(xi-X2) 



(22) 



(23) 



by virtue of the virial theorem (the power divergent self- 
energy graphs Fig.|2Ib), Fig.|2Ic), which renormalize the 
particle masses, vanish when evaluated in dimensional 
regularization). It can be shown in general that all other 
contributions to Sef / [xa] are down by powers of v relative 
to the particle kinetic terms or the Newton potential (see 
the appendix). It is then easy to deal with the presence 
of the large parameter L in the power counting. We must 
treat the operators that scale as Lv^ non-perturbatively, 
or in other words, we must expand about background 
fields Xq that to leading order in v satisfy the equations 
of motion of Newtonian gravity, treating operators which 
scale as iw" perturbatively. This is the resolution of the 
problem raised in |2,3| | regarding the apparent breakdown 
of a manifest perturbative expansion when the particle 
masses are larger than mpi. 

One can also show that diagrams that contain gravi- 
ton loops are supressed by powers of 1/L relative to the 
tree diagrams (after particle worldlines are stripped off). 
This is in accord with the intuition that loop diagrams 
represent quantum corrections which should not be kept 
in the perturbative expansion of classical field theory ob- 
servables. The fact that loop corrections are irrelevant 
to the quantities of interest also explains why we did not 
need to keep ghost terms in the functional integral of 
Eq. ini). 

What is the resulting theory after all potential gravi- 
tons are integrated out? It is a theory of radiation gravi- 
tons coupled to moments of the two-particle distribution 
in a gauge invariant manner. Working in the CM frame, 
J2a iTT-a^a = 0, the leading order contributions to pro- 
cesses in NRGR can be calculated from 
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FIG. 2: Diagrams contributing to Lu" potentials. The self-energy graphs in (b), (c 
effect. 



are pure counterterm and have no physical 
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where we suppress the argument {x'^, Xcm) of the canon- 
ically normalized radiation field h^,^, and write m — 

that after integrating out the potential modes, the only 
dependence on the orbital scale r appears in the Wil- 
son coefficients (the moments) of the theory. Thus in 
our theory, the decoupling between near-zone NR physics 
and far-zone radiation is manifest at the level of the 
NRGR Lagrangian. This explicit separation of scales, 
which is what motivates the EFT approach, should sim- 
plify the evaluation of higher order radiative tail effects 
within our formalism. Of the terms shown here, only the 
coupling of the (linearized) Riemann tensor to the mo- 
ment lij — ^^maXajXaj, Can source gravitational radi- 
ation, since all other terms in this Lagrangian couple to 
conserved quantities and therefore cannot radiate, up to 
higher order radiation reaction effects. In fact, the trace 
of the linearized Riemann tensor 
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mo 
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(25) 



2-,- ^dfh + OCh^) 

vanishes for on-shell graviton matrix elements, so that ra- 
diation only couples to the traceless quadrupole moment 
of the source 



Q^, 



(26) 



In the next section we will show more explicitly how to 
derive Eq. (|^ . 

Given Eq. (|24|l . how do we obtain the power radiated 
in gravitational waves from this formalism? We simply 
calculate the imaginary part of the self-energy diagram 
of Fig. 13 
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dx;fix^/y(a;;)/fei(x^) (27) 




FIG. 3: NRGR diagram whose imaginary part gives rise to 
the quadrupolar gravitational radiation power spectrum. 



X {RoiOj{Xi,^cm)RokOl (2^2' Xcm))- 

(We have only kept the part of the diagram that arises 
from the coupling lijRoioj, since it is the only term at 
this order which has a non-zero imaginary part). Com- 
puting the two-point function of the linearized -RoiOj one 
finds that it is proportional to the projection operator 
onto symmetric and traceless two-index spatial tensors, 
^^ikSji + ■^SiiSjk — ■^SijSki- Writing the above equation in 
terms of the Fourier transform of the quadrupole moment 



(fco) = / da:°Q,,(x")e-*'=°^°, 



we end up with 

Fig.El: 



(27r)4 fc2 + ie 



\QM\' 



(28) 
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This diagram then gives rise to the following contribu- 
tion to ImSef f[Xa] 
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80m% 
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.k4|Q,,(|k|)|2. (30) 



The integrand in this equation is proportional to the dif- 
ferential graviton emission rate over the history of the 
two-particle system. To get the differential power, multi- 
ply the integrand by an additional factor of the graviton 
energy |k|. In particular, including this extra factor of 
the energy and integrating over all momenta, we obtain 
the total power radiated in gravitational waves 
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(31) 



where the dots denote time derivatives, and the brackets 
denote a time average. Of course, this equation is just 
the usual quadrupole radiation formula. 



D. Examples at NLO 

The construction of the EFT in the previous section 
was perhaps somewhat formal, so it is helpful to see ex- 
amples of how to explicitly calculate the terms in the 
NRGR action. We will work out two textbook examples 
in this section: the Lagrangian that leads to the lead- 
ing order predictions for the power radiated in gravita- 
tional waves, Eq. (|24|) . and the EIH Lagrangian that 
describes the 0{v'^) general relativity corrections to the 
gravitational interaction of two point masses. 



1. EIH Lagrangian 

First consider Leih- It arises from the diagrams 
shown in Fig. 0| and Fig. |S1 Using the NRGR Feynman 
rules, the sum of the diagrams with a single potential 
graviton exchange are given by 
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with X12 = Xi — X2 
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and 



Fig. lit = [dhdt2Sih-t2 



-lk-(xi-X2) 



3i 



GNmim2 2 
dt- rvf. 



2 J |Xi - X2 

Here, we have used the formula 



(34) 



1 T{d/2-a) fx 



(35) 

However, at this order in u, the single exchange dia- 
grams are not sufficient. The NRGR power counting tells 
us that we also get contributions from diagrams with two 
and three graviton propagators. For example, the graph 
in Fig. ^ scales as 
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The first factor stems from three insertions of leading Einstein-Hilbert Lagrangian to order H^, it is straight- 
order vertices while the second bracket arises from an forward to show 
insertion of a three-graviton interaction. Expanding the 



(T«(:r0i?£°(x2)7?^°(x3))> =-£^^(x;-x°)^(.;-x^)(27r)3<53 ) IT ]^ x E (37) 

\ r / r— 1 ^ r— 1 

SO that 

Fig.EK = f / e^i:.''--(r«(^i)i?^^(^2X(a;20)) - - J rfi^r^vlv^' (38) 

where we have dropped linearly UV divergent self-energy contributions which can be absorbed into the particle masses, 
and which vanish when evaluated in dimensional regularization. Finally, the seagull diagram, Fig.jSJ), also contributes 
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FIG. 4: Diagrams contributing to Lv'^ terms in the two-body potential. The ® in (a) denotes an insertion of the potential 
graviton kinetic term. A similar diagram to (c) with 1 +-> 2 is not shown. 




(a) (b) 

FIG. 5: Diagrams with two and three potential graviton lines that also contribute to the Lu^potentials. Similar diagrams with 
1^2 not shown. 



at order Lv^, and is given by 
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Putting together Eq. (|^ . Eq. Eq. ||23J, Eq. and Eq. together with their mirror images under the 

interchange of the particle labels, we end up with 
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(40) 



where we have also included the relativistic corrections to 
the kinetic energy of the point particles, which contribute 
to the effective action at order Lv^. This Lagrangian 
was first derived by Einstein, Infeld and Hoffmann 
by different methods. A diagrammatic derivation of the 
static part of this Lagrangian, using Lorentz co- 

variant worldline methods [2^ . was given in ref. [29| . A 
quantum field theoretic derivation of the 1 /r^ potential, 
treating the source particles as dynamical fields, was pre- 
sented in m 113. 



2. Couplings of Radiation Gravitons 

It is also instructive to use the formalism to derive 
the leading order Lagrangian that describes the interac- 
tions of radiation gravitons with NR sources, Eq (|24|l . 
The derivation involves all aspects of the construction 
of NRGR (the multipole expansion for radiation, inte- 
grating out potential gravitons, velocity power counting 



rules). At the lowest two orders in the expansion, the in- 
teraction of radiation with the NR sources follows simply 
from the multipole expansion of the radiation field. Ex- 
panding the proper time time term in powers of u, using 
the counting rules in Table I, leads to 

Lyi/2 = -^r^y^mahoo. (41) 
2mpi ^ 

a 

L^3/2 —Pc^iJlQ, = 0, (42) 

mpi 

in the CM frame. Note that to preserve manifest scal- 
ing in V we must multipole expand about the center of 
mass, as defined below Eq. (|24|l . These terms (when in- 
tegrated over time to form an action) scale as ^/Lv^^"^, 
and \fLv'^l'^ respectively according to the power count- 
ing rules of the EFT, and generate diagrams which (i.e. 
terms in Sf^ff\XcS) that are proportional to a single power 
of i. 

At the next order in the velocity, the situation is not as 
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straightforward. Keeping terms with one more power of 
V (either from the multipole expansion or from expUcit 
powers of the particle velocities) one finds, in the CM 
frame, 



-'mult 
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This Lagrangian has two problems. First, it is not gauge 
invariant under infinitesimal coordinate transformations. 
Second, it seems to predict that unphysical 00, Oi gravi- 
ton polarizations can be emitted by the two-particle 



bound state. Of course, these two problems are related. 
The resolution to this is that at this order in the expan- 
sion there are additional contributions to the effective 
Lagrangian, from diagrams with both external radiation 
gravitons and internal potentials, such as those of Fig.|Sl 
The diagram in Fig. 03 is given by 



Fig.EjD 



27 ^ |X1-X2| ' 



(44) 



For the diagram in Fig. IHt, we need the hH^, vertex. It 
can be obtained from 
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where the integral over momentum k has been sup- 
pressed, and because of the multipole expansion, h^^ is 
evaluated at (a;°,Xcm)- Given this term, it is easy to 
show that 
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as well as the equations of motion for xi^2(i) from the 
order Lv° effective action, we find 
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Adding together the result of Eq. (|44(l plus its mirror 
image under 1^2 and Eq. to Eq. we find the 
complete NRGR action at order 
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There are several things to notice about this equation. 
First, we find that /iqo couples at this order in v to the 



Newtonian gravitational energy of the bound system. At 
this order in the expansion this is a constant of the mo- 
tion, as is the orbital angular momentum, which couples 
to hiQ. Consequently, neither of these quantities can ra- 
diate. This is in accord with the expectation that hoo,hio 
do not correspond to physical graviton polarizations, and 
should not be emitted by a graviton source. It should be 
possible to give a general diagrammatic proof that the 
diagrams in the matching involving one external Hqq are 
related in a simple way to the diagrams with no external 
radiation gravitons, i.e, the diagrams that contribute to 
the total energy, as required by general covariance. 

Also, note that this Lagrangian is manisfestly gauge in- 
variant under infinitesimal gauge transformations. This 
is obvious for the last term in Eq. (|49|l . which involves 
the (linearized) Ricmann tensor. For the first term 
use the fact that the Newtonian energy is (at this or- 
der) time independent, so that under the transforma- 
tion hoo — > hoo + 2do^o, this term in the Lagrangian 
transforms into a total time derivative. Likewise, under 
hio —>■ hio + di^o + do£,i, the coupling of hio to the angular 
momentum transforms into a total time derivative. We 
expect that the requirement of gauge invariance should 
provide strong constraints on the form of the NRGR La- 
grangian. In particular, the derivation of the 0{v'^) ra- 
diation couplings presented here could have been signif- 
icantly shortened had we imposed gauge invariance. In 
practice, one can reduce the amount of work by keeping 
only the terms from the multipole expansion and fixing 
the remaining terms by demanding gauge invariance of 
the radiation couplings. 

It is interesting to note that in order for the couplings 
of the radiation gravitons to preserve gauge invariance. 



1 




T 



I 
I 

2 , . i ► , i . 

(a) (b) (c) 

FIG. 6: Diagrams that contribute to the matching of radiation graviton couplings. The diagram (a) is the multipole expansion 
carried out to order w^. Diagrams (b), and (c) are needed to ensure gauge invariance of the radiation graviton couphngs. 
Similar diagrams with 1 <-> 2 not shown. 



it is crucial that we include the diagrams of Fig. |SJd and 
Fig. El;. These diagrams, which arise due to the non- 
Abelian nature of the gravitational interactions, repre- 
sent the fact that in gravity it is not just the ensemble of 
NR particles that acts as a source of gravitons, but also 
the energy stored in the gravitational field that it sets 
up. To obtain consistent results such non-linear effects 
must be included. Without the velocity power counting 
rules of the EFT, it would have been difficult to ascertain 
whether such non-linear diagrams contribute already at 
this order in the velocity expansion. In the classical gen- 
eral relativity approach to calculating the power emitted 
in gravitational waves by a NR source, this non-linearity 
appears in the expression for the far zone gravitational 
field as an expansion in moments not over the stress ten- 
sor T^'^ of the matter alone, but also over a suitable (non- 
unique) stress-energy "pseudotensor" r^" for the gravi- 
tational field itself. This pseudotensor is defined in such 
a way that d^T^^ + d^T^^ = 0, and physically plays the 
same role as the diagrams of Fig. EId and Fig. ^ do in 
our formalism. As an aside, we also note that the ne- 
cessity to include diagrams involving gauge particle self- 
interactions in order to derive the couplings of radiation 
to sources has an analog in the worldline approach to 
NR bound states in QCD [sj. For instance, to show 
that the color dipole moment of a QQ bound state cou- 
ples to the full color electric field, it is not enough to keep 
diagrams with gluons emitted from the quark worldlines. 
One also needs to keep diagrams involving the 3-gluon 
vertex, which can be shown to contribute at the same 
order in the velocity [s^ . 



III. DIVERGENCES, RG FLOWS AND FINITE 
SIZE EFFECTS 



The coefficients of the relativistic point particle theory 
of Eq. (jSJ exhibit non-trivial RG flows, even at the clas- 
sical level (RG flows in classical fleld theories coupled to 
extended spatial sources have been analyzed in |33 ) . The 
non-trivial RG structure of the theory is associated with 
logarithmic UV divergences which arise due to the singu- 
lar nature of the point particle limit. To renormalize the 



theory we consider a single particle at rest at x = and 
compute the diagrams that contribute to the terms in the 
effective action for a background gravitational field that 
are linear in the background /i^^ = — 77^^ . If we ignore 
for now the non-minimal curvature couplings of Eq. jSJl, 
this is given by diagrams such as those of Fig. [7| with a 
single external graviton and any number of insertions of 
the worldline proper time coupling, each proportional to 
the particle mass (we restrict ourselves to tree diagrams, 
which are the ones that remain in the limit mpir 3> 1, 
with r = |x|). 

While the diagram with a single mass insertion is fi- 
nite, diagrams with more mass insertions have UV diver- 
gences. For instance the diagram with two powers of m in 
Fig.[7fb) contains a piece which has a power divergence of 
the form ~ J d'^~^q/q^ , as can be seen by counting pow- 
ers of momentum from the internal propagators and the 
graviton self-interaction vertex. An explicit calculation 
shows that this linearly divergent piece, which vanishes 
in dimensional regularization, just renormalizes the bare 
mass parameter m appearing in the point particle ac- 
tion. Besides this power divergent piece, the diagram of 
Fig. Cfb) also gives rise to a finite, non-analytic term in 
the background field effective action "sS] of the form 



(50) 



where, in background field gauge with a gauge fixing term 
for the "quantum" graviton field H^^ of the form Sg f = 
/ d^x^r^F'S F^ = - \D^'H'^ one finds 



T(2)^.(fc) = (27r)^(fc") 
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(51) 



where = (1,0) is the particle four-velocity. As ex- 
pected by background field gauge invariance, we have 
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FIG. 7: Some contributions to the background graviton effective action in the relativistic point particle theory. Diagram (b) 
has a power divergent term. Diagrams (c) and (d) are both power and logarithmically divergent. Two other diagrams that 
arise as permutations of (c) are not shown. 



kfjT^^^{k) — 0. In fact, since r[(7^^] must be covariant 
with respect to the background gravitational field, the 
above result is sufficient to completely fix all terms in 
the effective action that are proportional to two powers 
of m and any number of powers of the background field 

Although the power divergences that arise in the cal- 
culation of T^^{k) have no physical consequences, log- 
arithmic divergences can give rise to measurable effects 
since by dimensional analysis, in a massless field theory 



the logarithm of the cutoff is always accompanied by a 
logarithm of the external momentum. Consider the di- 
agrams with three insertions of the particle mass. They 
give rise to a term in the effective action of the form 



ng..] = ^J ^,h,A~k)Ti-;ik) + ■■■, (52) 



where we find, from Fig. Oc) and Fig. [Tfd) 
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As a non-trivial check of this result, we note that 
gauge invariance with respect to the background field, 
A;^r^3^(fc) — 0, does not hold diagram by diagram, but is 
only manifest after the terms in Fig. [7|[c) and Fig. Od) 
are combined. To obtain this equation, we worked out the 
necessary background field gauge Feynman rules [2^ in d 
dimensions, and performed an expansion in e = {A — d)/2, 
keeping only the terms that do not vanish as e — s- 0. 
Tf^^^ik) contains a UV logarithmic divergence, as can be 
seen by power counting the integral 

°^ ' J (27r)'i-i (27r)'i-i q2p2(q + p + k)2' 
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which evaluates to 
/o(k) = 
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(55) 



In the second line we have expanded about d = 4— 2e and 
introduced an arbitrary subtraction scale /i. Inserting 



this expression back into Eq. H53|l . we see that the set of 
diagram with three mass insertions has a 1/e pole which 
cannot be absorbed by a shift in the mass parameter m 
of the minimal point particle action, whose Feynman rule 
is, for a particle moving on a straight line trajectory with 
constant four- velocity , 
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contribute uchhd uu _l i^y^^j 
{1/mpi) J d^k/i2nyh,A^k)Tl^-{k), 
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of Eq. © 

the form 



with 



T^7(/c) = (27r)(5(fcO) CR{T]'"'k'^ - k^'k") + ^cyk^v^'v 

(57) 

where we have dropped terms proportional to fc" which 
vanish due to the delta function. Because the divergent 
terms in Eq. 1)53(1 have the same dependence on k and v 
as the contribution from the curvature couplings, we see 
that the 1/e poles can be subtracted away by a proper 
choice of coefficients cpy. In the MS scheme, the renor- 
malized 0(rn^h) term in the effective action now becomes 
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(58) 



where CR^vifJ-) are the renormahzed couphngs. 

Since the renormahzation scale /i introduced above is 
arbitrary, we must have fidT[g^i,]/dfi = (since the clas- 
sical field /i^^ does not pick up an anomalous dimension). 
Thus the explicit dependence of the logarithm in Eq. H58() 
on the subtraction scale fj, must be cancelled by allowing 
the coefficients cr^v to vary with scale. Our theory there- 
fore exhibits non-trivial classical RG scaling, 



dcR 
dfi 
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As a consistency check, we may use the linear term in 
r[g;^j/] to calculate the spacetime metric generated by a 
point mass at rest, up to order (GNm/r)^ . Fixing the 
gauge r/^T^^ 
simply given by .g^,^ 
space 



for the background field /i^^, this is 



h^y, where. 
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h=0 



Using our results for Tl^^^{k)^ we find, dropping delta 
function contributions which first arise at order m? , 
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and goi = 0. It is straightforward to show that this metric 
satistfies the vacuum Einstein equations = away 
from the point x = where the point particle source 
is located, independently of the value of the renormal- 
ized parameter c/j(/Lt). This is consistent with Birkhoff's 
theorem, which states that the space of static, rotation 
invariant solutions of the vacuum Einstein equations is 
one-dimensional, parametrized by the mass m. Alterna- 
tively, one may check the validity of our result for g^^ by 
constructing a coordinate transformation to the more fa- 
miliar harmonic coordinates x'', with 5^'^(a;)r^^(a;) — 0, 
which in terms of our coordinates x'^ are given by, up to 
0{Gjfm'^ /r"^) terms. 



(64) 



i+2[^y-u^y (65) 
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Thus, as required by general arguments, the parameters 
Cfi vifJ-) a'l'e not observable in the long range field of iso- 
lated static sources. 



As we argue in the next section, it is likely that the 
value of the matching coefficients cjiydio) are themselves 
of order Gj^m^. In this case, it is easy to see using the 
NRGR power counting rules that diagrams with one in- 
sertion of the non- minimal operators in Eq. jSJl, for in- 
stance the diagram of Fig.lSJa) lead to corrections to the 
two-body interaction Lagrangian which are suppressed 
by a power of relative to the leading Newton interac- 
tion. However, it can be shown that the worldline op- 
erators J drR, J drR^iyX^x'^ can be removed by a local 
redefinition of the metric tensor (that is, the coefficients 
c_R,v can be shifted arbitrarily by these transformations). 
Given that any possible signature of internal structure 
that may arise at must be encoded in these two op- 
erators withing a point particle description, we conclude 
for spinless objects there are no finite size effects up to 
order . 

Even though all values of fi give rise to the same 
physics, in practice certain choices may be more conve- 
nient when doing actual computations. In general, this 
leads to renormahzation group equations which can be 
used to simplify the calculation of logarithmically en- 
hanced terms in the perturbative expansion. We will 
give an example in Section of how to exploit the free- 
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dom in choosing /i for the unphysical couplings cv,r to 
simphfy the calculation of certain contributions to the 
two-body potentials that are enhanced by a logarithm 
of a = rs/r, where is a length scale that character- 
izes the size of the extended object. The calculation of 
logarithms for physical tidal couplings, which involve the 
Riemann tensor, is completely analogous to this case. 

Any choice for the couplings Cfl_y(/i) can be trivially 
expressed in terms of the parameters at some fixed scale 
Ho by solving the above RG equations 

6 Ho 
cvijj) = cv{ho) + i^G^iTir'hi—. (66) 

6 Ho 

(as we will see later, at higher order in the velocity expan- 
sion, the coupling constant flows show more interesting 
behavior than is depicted in these equations). If the cou- 
plings CR^vifJ-o) had any physical content, in order to fix 
them say at a scale of order the characteristic size Ts, 
requires model-dependent information about the details 
of the short wavelength physics that resolves the point 
particle divergences. This will be true for the physical 
tidal couplings as well. This physics is independent of 
the NR two-body dynamics, so that no PN formalism 
based on point particles can, by itself, give rise to a full 
prediction for the gravitational wave observables. In the 
next section we will explain how, given this short distance 
model, it is possible to perform a matching calculation 
that fixes the values of the worldline operator coefficients 
in terms of the parameters (size, equation of state, . . .) 
that characterize the extended object. 

An interesting question arises as to the size of the 
matching coefficient cr,v{ho), or their more physical 
tidal counterparts. If it turned out that these coefficients 
were anomalously large^, then this would put the system- 
atics of the PN approximation in jeopardy. For instance 
if the coefficients in Eq. ||SJ) were much larger than Gj^m^ 
(no explicit powers of v, a dynamical variable, can appear 
in operator coefficients), the Newtonian result would no 
longer be leading order, and the PN expansion would 
cease to be valid. Although it seems rather unlikely that 
this constraint is violated, it is interesting to note that 
within the EFT formalism here, we would retain calcula- 
bility given a model for the short distance physics, but we 
would have to change our power counting scheme in order 
to accomodate the lack of a Newtonian gravity limit. 



IV. MATCHING ON TO THE POINT 
PARTICLE THEORY 

A. Outline of the Matching Procedure 

If one works at high enough of an order for the world- 
line operators to become relevant, i.e. beyond v^, then 
given that PN theory has nothing to say about the actual 
values of the worldline operator coefficients, one needs a 
model of the internal structure. Otherwise the best one 
can do is to express the predictions for the gravitational 
wave observables in terms of the parameters of the EFT 
and use the gravitational wave data to obtain information 
about the magnitude of the coefficients. Alternatively, 
if one does have information about the structure of the 
binary constituents, one can fix the values of the cou- 
plings in the effective action at a renormalization point 
fi ^ l/rg. In this case the data can, in principle, be used 
to test one's assumptions regarding the structure of the 
binary star sources. In this paper we will not perform the 
matching calculation, but instead just sketch how such 
a calculation would lead to an unambiguous prediction 
for the gravitational wave observables. We will illustrate 
how this works within the context of a specific toy model 
which captures the essential physics. 

The matching procedure is standard within effective 
theories , but to keep the discussion self-contained we 
will review here the basic concepts. The idea is simply 
to compare observables such as Green's functions in the 
"full" (valid at short distances) and effective theories and 
adjust the value of the EFT couplings in such a way 
that the effective theory reproduces the results of the 
full theory. The effective theory is designed to reproduce 
the long distance physics, characterized by non-analytic 
behavior about zero external momentum, but it may get 
the short distance, which is analytic in the momentum 
(that is, local in coordinate space), wrong (hence the 
divergences) . 

Consider for instance the matching of one-point func- 
tions in the full theory (consisting of an extended source 
interacting with gravitons) and the EFT. In the full the- 
ory, the one point-function of the graviton are of the form 

G^'\k^)^^F{k\K,), (67) 

where F{k'^,Ki) is in general a complicated function of 
the external momentum k as well as some dimensionful 
parameters which describe internal structure. For sim- 
plicity we assume in the discussion that the only relevant 
scale is the size of the object . The EFT is valid in the 
kinematic regime k^r^ ^ 1, so in order to match, one 
expands the function F{k'^r'j.) about this point. In the 
simplest situations, this expansion takes the form 

G^'\k') = i [P{k\l)+Q{erl)Herl)] , (68) 



An explicit example of this arises in some mo dels of branes em- 
bedded in higher-dimensional spacetimes l35| . 



where P{z),Q{z) are polynomials. In general, the non- 
analytic dependence on the momentum may be more 
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complicated, involving polynomials times more general 
functions of In/c^r^, however, we focus on this case for 
simplicity. In practice, since we are interested in a low 
energy expansion, we need to keep only the first few terms 
in the polynomials P{z),Q{z). Provided that the point 
particle coupled to gravity properly reproduces the IR 
(long distance) physics, the EFT result must be of the 
same general form as the above 



p{e,c^{^l)) + Q{k' 



.Cj(M))ln — 



(69) 

where the parameters Ci (fj.) are the coupling constants of 
the EFT defined by some suitable renormalization pre- 
scription, which in the point particle theory are the parti- 
cle mass and the curvature couplings. The functions P, Q 
are polynomials of k^ whose coefficients are products of 
the Ci (/i) . Note that the log divergence arises because we 
have taken the singular limit 0. Since this loga- 

rithm arises in the full theory from the light degrees of 
freedom which the EFT must also contain, it is crucial 
that the non-analytic momentum dependence in Eq. H69|l 
coincides with that of Eq. (|^ . In order for the EFT to 
reproduce the results of the full theory (at least to some 
fixed order in the small parameter A;^r^), we then tune 
the coefficients cji^vip-) so that G^^-'(fc^) — G^]^]pj,{k^) in 



EFT\ 

we had P ^ 
'/mpi, hence 



the hmit fc^r^ < 1. 

For instance, in the previous section, 
m/mpi + CRy{fi)k'^ /mpi, and Q ^ G%m^ 
in the full theory that resolves the point particle singu- 
larities, we expect Q{z) ~ G%m^z/r'^mpi for z ^ (we 
are ignoring finite, non-analytic power corrections of the 
form \/k^ from FigEbi which are not relevant to this dis- 
cussion). Comparing Eq. (|58f) and Eq. H69|) in this case 
we obtain a relation 



CR,v{p) ^ Co + G%m^ In(^rs), 



(70) 



where cq is a full theory parameter obtained from the ex- 
pansion of P{z) to linear order. To avoid large logs which 
could hinder our perturbative expansion, we choose /i ~ 
l/rs- This is the boundary value of the solution to the 
RG equation Eq. H66|l . Actually, this equation ignores the 
presence of scheme dependent constants (which are not 
relevant to physical predictions) which are of the same 
magnitude as the coefficient of the logarithm. Thus even 
if the short distance parameter cq were much smaller than 



G 



N' 



we would find cpy{fj,o) ~ G 



N' 



B. A Toy Model of the Short Distance Physics 

In this section we illustrate the general remarks above 
in the context of two simple toy models. The matching of 
the worldline operator coefficients which carry informa- 
tion about the static properties is conceptually similar to 
the matching in a model consisting of a self-interacting 
scalar field coupled to a point-particle source. To match 



point particle couplings which carry dynamical informa- 
tion about the binary star constituents, for instance a 
tidal operator of the form J drRfj.^apR'^'^'^^ that may 
appear in the point particle EFT, requires matching am- 
plitudes with more than one external graviton. We will 
consider an analog of this in electrodynamics, where we 
will describe the matching onto the point particle opera- 
tors / d.TFf.^Ft"', J dTxt'x''Ff,aF°"' which encode the dy- 
namical repsonse of a charge distribution to applied low 
frequency electromagnetic fields. Explicit calculations of 
the counterterms need to fix the operators coefficients for 
the point particle EFT a gravitating compact object are 
currently in progress |l9l |. 

First consider a self-interacting massless scalar field 
4>{x) in four dimensions coupled to a localized source 

C^l{dcj^)'-^cb'-pix)q^ix), (71) 

and imagine that through some unspecified dynamics, 
there is a frame in which the source p{x) acquires a profile 
of the form 



p(x",x) 



9 



^3/2„3 



exp 



(72) 



where and Vs are constants that parametrize the short 
distance physics. Although this model is unrealistic in 
the sense that the profile p{x) does not arise dynami- 
cally but is rather put in by hand, it is adequate for our 
purposes here. 

In the limit where only long wavelength modes of 
the field (j){x) are accessible, the detailed structure of 
the source p{x) becomes irrelevant. As far as long dis- 
tance modes are concerned, one may replace the localized 
source with a point particle that has suitable local cou- 
plings to the field (j){x). 



Seff — S[(t)] — TO J dr ~ g J drcp 



(73) 



where S[(j)\ is the bulk scalar field action, and we have 
suppressed higher dimension worldline couplings involv- 
ing more powers of (j) or its derivatives. The coupling 
g appearing in this equation is a function of the para- 
maters g, rg of the short distance source model, which 
we determine by matching amplitudes in the low energy 
EFT with those obtained from Eq. (|71|) . At lowest order 
in A, this is trivial. Consider for instance the amplitude 
for vacuum —t (f>. In the full theory 



i{2TT)S{k°)Afuii = -t / d^xp{x){k\cj){x)\0) (74) 



so that 



(75) 



which essentially is the form factor associated with the 
source distribution p{x). In the EFT, 



t{2Tr)S{k°)A 



EFT 



-ig J dT{k\(t){x)\0) 
^z{2n)S{k")g. (76) 
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The two amphtudes agree in the Hmit k^r^ if we 
set g = g. At higher order in the momentum expansion 
we would also need operators of the form J drd'^^'cj) in 
the low energy theory to reproduce the full theory am- 
plitude. Similarly, by matching the vacuum to vacuum 
amplitudes, we find m — (l/4\/27r)5^/7's- 

Note that in the point particle EFT, there are logarith- 
mic divergences in the vacuum —> (j) amplitude at 0{g^X), 
from a diagram analogous to the one in Fig.Ud), 



iA 



EFT 



~ig 



iXg^ 
3! 



(77) 



-«5(m) 



1927r2 



In — 



where /o(k) was defined in the previous section, and as 
before we have used dimensional regularization plus the 
MS scheme to renormalize the divergent amplitude. If 
the EFT is to be a correct description of the low energy 
physics, the non-analytic dependence on the momentum 
in this equation should reproduce a similar log in the full 
theory amplitude. Indeed, calculating in the full theory 
and including insertions of the (t>* interaction we find, 
dropping terms subleading in k^ 



0^ 
r? < 1 



iA 



full 



-ig- 



-ig 



iXg^ 



1927r2 



(2^)3 q2|q 



2J2 



In- 



k^r 



k| 

7-1 



(78) 



which agrees with the claimed universal result Eq. (|68|l . 
Therefore 



1927r2 



In 4^+7 + 2 



(79) 



which indicates that to avoid large logarithms that may 
render perturbation theory invalid, we should match at 
a scale fi ~ ^/ts, so that g{l/rs) ~ g. Then we can use 
the RG equation 



dg_ 



967r2 



(80) 



to sum logarithms of k^r^. Even though we phrased the 
matching problem using quantum field theory language, 
we stress that our discussion here applies also in the clas- 
sical regime. This can be seen by matching not scattering 
amphtudes, but the 1-point functions {(j){x)), in the pres- 
ence of the point particle source, which have an obvious 
classical interpretation. 

How does this model bear on the matching problem 
for the point particle coupled to gravity? As in the grav- 
itational case, our toy scalar field model exhibits UV di- 
vergences in the point particle limit. In both cases, we 
expect those divergences to be resolved by the finite ex- 
tent of the distribution which sources the field. We saw 
this explicitly in the scalar field theory, where we learned 
that the matching procedure fixes the couplings of the 



point particle at a length scale of order the size of the 
source. We expect that the same will occur in the gravi- 
tational problem. 

However, there are several aspects of the gravity prob- 
lem that are not captured by this simple model. In the 
scalar field theory, the UV divergences associated with 
the point particle singularities affect the leading order 
coupling J drcj) itself, and thus the point g = is a (triv- 
ial) fixed point in the RG sense. In the gravitational case, 
on the other hand, the structure of the RG flow is such 
that even if the curvature couplings were equal to zero 
at /J, ~ l/fs, non-zero values would get induced at longer 
wavelengths. Second, in the scalar model, the order A^^ 
correction to the long range field encodes information 
about the internal structure of the extended object that 
couples to (j). For a source in isolation coupled to gravity, 
Birkhoff's theorem states that the curvature couplings 
cannot generate long range static forces that could carry 
information about the short distance dynamics. Never- 
theless, Birkhoff's theorem does not prohibit the momen- 
tum space off-shell graviton one-point function from hav- 
ing sensitivity to the internal structure. Thus in principle 
one could still match the graviton one-point functions (or 
the quantities T^^[k) appearing in the background field 

action), although in practice it may be difficult to com- 
pare in a gauge invariant manner. Because of this, it 
may be simpler to match by comparing gauge invariant 
quantities, such as a the amplitudes for fields to scat- 
ter off the geometry of the extended source versus the 
scattering amplitude in the point particle theory. 

Fixing the coefficients of the operators that measure 
tidal effects in the point particle EFT will require match- 
ing amplitudes with more external gravitons. A simi- 
lar situation arises in electrodynamics, where dynamical 
effects such as induced dipole moments are encoded in 
the photon two-point function (equivalently the response 
to an external field). Consider for instance the effec- 
tive point particle action that describes the interactions 
of the electromagnetic field with a perfectly conducting 
grounded sphere of radius Ts- Since by definition the 
times scales associated with the dynamics of a perfect 
conductor are short compared to those of the electromag- 
netic fields it interacts with, the electrodynamic analog 
of Eq. lO is then 



dr + eQ / dxf'A 



(81) 



Imagine putting this system in a constant background 
electromagnetic field T^j,^. Expanding the field about T 
in the point particle EFT, the fiuctuating photon field 
A^ acquires a vacuum expectation value which is given, 
in the Feynman gauge, by 



2^b^^''^^A^'^^'W^^''- ^^^^ 
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FIG. 8: Sample contribution to two-body potentials at In a. 
All particle vertices are the leading order mass «° mass inser- 
tion. The ® denotes an insertion of the potential graviton 
kinetic term. 

Imposing the appropriate boundary conditions for a per- 
fect conductor, the "microscopic" model predicts 

so matching fixes the gauge invariant EFT coefficients to 
be Q = 0, /3 = 6a = 67rrf . In this particular model, all 
higher dimension operators built out of more derivatives 
or powers of have vanishing coefficients. 

V. LOGARITHMICALLY ENHANCED 
POTENTIALS 

Although the operators J drR and J drx'^i'^ R^^i, can 
be removed by a redefinition of the metric and thus never 
contribute to physical observables, it is instructive to use 
these operators to illustrate the general procedure for 
calculating the contributions to the two-body potential 
which are enhanced by a logarithm of a = Tsjr. Such 
terms are independent of the details of the short distace 
dynamics. Although these terms have been obtained in a 
different gauge by conventional PN methods in Refs. , 
here we give a derivation that exploits the RG scaling of 
the coefficients in the EFT to simplify the calculations. 

The idea is that because the amplitudes in the EFT 
are independent of the subtraction scale /i we can choose 
any convenient value for it. In particular, if one chooses 



/i to be of order the orbital scale l/r, then the terms in 
the two-body potential that involve logarithms of /i|xi2| 
(which arise from taking the Fourier transform of dia- 
grams involving subgraphs such as those in Fig. Efc,d) 
are minimized. In this case the log enhanced contri- 
bution to any amplitude can be obtained entirely from 
diagrams with insertions of the local operators J dri?, 
/ drxf^i'^ Rfj^^, with coefficients cuy evaluated at a scale 
fi ~ l/|xi2|. This procedure amounts essentially to re- 
moving the short distance modes of the gravitational field 
with wavelengths between the internal distance rg and 
the orbital scale r. 

For example, consider the diagram shown in Fig. |S1 
which contributes to the In a terms in the potential. 
Schematically, it contains terms of the form 

Fig.El^ J dx°-^^^^[c{^^) + G%mlH^^\^U] (84) 

where c{fi) denotes the contributions that arise due to in- 
sertions of the non-minimal curvature couplings and the 
logarithmic piece arises from a calculation similar to that 
of the 0{m^) terms of the background graviton effective 
action in Section Uni Altho ugh it is hard in general to 
calculate those logarithmic pieces, their argument is in- 
variably of the form /z|xi2|, so by choosing fi ~ l/|xi2| 
we may eliminate these pieces in favor of the effective 
couplings c(l/|xi2|). These in turn may be obtained in 
terms of the bare parameters c(l/rs) by using the RG 
equations for the couplings in cj^y in Eq. (|59|l . Eq. H60|l . 

c(l/r) - c(l/r,,) + G^m3 In a. (85) 

Since the validity of the velocity expansion requires 
c(l/rs) < G^m'^, the logarithmic term gives the domi- 
nant contribution to the two-body potential at this order 
in V. 

Following this line of reasoning, we can obtain the 
model independent In a terms by calculating the di- 
agrams of Fig. |51 Fig. ^1 and evaluating the coefficients 
CR.vifJ^) at the scale fi ^ l/|xi2|. We find 



In a 



AGn (m24'^(|xi2|) + mi4'^(|xi2| 
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4GArm24^(|xi2|) 



X12 
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3(V2 ■X12)' 
|X12P 



Xl2 



2) 



4G^mim2 



X12 



(c«(|xi2|) + 4'^(|xi2|)) (86) 



In this equation, the first and second terms follow from 
the diagrams with a single graviton exchange in Fig. |51 
while the last term comes from the diagrams with two 
and three graviton propagators. Fig. ^| For the dia- 
grams involving insertions of the operator J drR^ijX^x'^ 



one finds that the one-graviton exchange diagrams van- 
ish identically while the diagrams from Fig. IIUI cancel, so 
that the two-body potential is independent of the coeffi- 
cients Cy . We have dropped parts of diagrams that give 
rise to a contact delta function potential that does not 



(a) 



(b) 



(c) 



FIG. 9: One potential graviton exchange diagrams contributing to the two-body Lagrangian at v^' \na. The heavy dot indicates 
an insertion of the J drR or J dTRn^x'^x" vertex with c_b,i/(a') evaluated at n ^ 1/r. Similar diagrams with 1 2 not shown. 



generate long range forces, as well as terms that can be 
expressed as a total time derivative and therefore do not 
contribute to the equations of motion. Note in particular 
that had we not chosen the subtraction scale /i to cancel 
the logarithms coming from the divergent integrals we 
would have needed, in order to obtain Eq. H86(l . a con- 
tribution from a diagram involving the quintic vertex in 
Hfj_i, from the Einstein-Hilbert Lagrangian. By choos- 
ing fj. ~ l/|xi2|, this piece of the potential arises instead 
from the diagram in Fig. II OF . which is much easier to 
calculate^. To see this note that the insertion of the part 
of the operator J drR that is quadratic in -ffk/^i/ in the 
diagram of Fig. [TUf is related by gauge invariance to the 
short distance piece (the divergent log) that is generated 
in Fig. □ 

Defining the binding energy of the binary system in 
terms of the two-body Lagrangian a,s E = J^a ' 
{dL/dVa) — L we find for a circular orbit of radius r 
in the CM frame m-aXa = 0, 



the metric in terms of a new metric g^^i, as 



1 ( Gi^m 



In a 



(87) 



where [i = m\milm^ m ~ mi + m2, and rg is a length 
scale of order the radii of the binary star constituents, 
which we assume for simplicity to be of similar size. 
One can check that this result agrees with the energy 
of a test particle of mass fi in circular motion around a 
Schwarzschild black hole of mass m (in the coordinates 
used in Section llll() in the limit — > 0. This requires 
knowledge of the logarithmic term in goo{r) proportional 

to l/r"*, goo = 1 H 1- |(GArTO/r)'* In fir. 

Because r is not a physical quantity, it is more natural 
to express the binding energy E in terms of the orbital 
frequency lu, which can be measured by observers at in- 
finity. A short calculation shows that the contribution 
to E{iLj) due to the coefficients CRyifj) exactly cancels, 
which by RG invariance also implies that the logarithmic 
dependence on /i must disappear. This follows from the 
metric field redefinitions alluded to earlier. If we write 



o 2 



cIt 



V9 



we find that to linear order in ^ 



-2m% 



d^x^R[g] = -2t 



'^Pi 



d'^x^/§R[g]+^ J drR, 
(89) 

so that the theory with curvature coupling cr is equiv- 
alent to a theory with coupling -t- ^. Since ^ can be 
chosen arbitrarily without affecting physical predictions, 
it can be set to ^ = —cr in which case the effects of the 
operator J drR are completely removable^. Similar field 
redefinitions can be performed to remove cy. This con- 
clusion can also be reached at the level of Eq. H86() . Since 
to derive Eq. (|86|l we had to fix a choice of gauge, a shift 
in the metric becomes equivalent to a gauge transforma- 
tion. It is then easy to see that if we let Xq ^ Xq + Sxa, 
(a — 1,2) with (5xi = — 4GArc^-'(/i)xi2/|xi2p and 6x2 = 

4GArc^''(^)xi2/|xi2p, the shift in the leading order New- 
tonian Lagrangian cancels the potentials calculated 
above. 

Beyond these spurious operators, any other operator 
involving the Ricci tensor is also removable by suitable 
metric redefinitions. This means that only operators 
built from the Weyl tensor have any physical content. 
The simplest physical operators that can be written down 
in a theory of spinless particles are then 

/ drG^.a^jG^''"/',/ rfri^i^G^.^pG^""- , 

/ dTx^xf^XpX^Cf.^.pCP"", (90) 

where we did not build any operators using the four- 
acceleration a^, since such terms can be set to zero by a 
redefinition of the point particle worldline. The NRGR 
power counting rules then indicate that these operators 
do not give finite size contributions to the gravitational 



This is true at least when the RG equations Eq. 1591 , Eq. 1601 are 
calculated in background field gauge, as in the previous section. 



Note however, that the shift in Eq. I88i induces in the two-body 
sector contact interactions of the form J drS^ {xi {t) — X2 (t)) /^/g 
which are physical, although irrelevant in classical processes. 



(d) (e) (f) (g) 

FIG. 10: Two and three potential graviton exchange diagrams contributing to the two-body Lagrangian at v^lna. Similar 
diagrams with 1^2 not shown. 



wave observables until at least order v^'^. This is consis- 
tent with Newtonian expectations. A complete account 
of finite size effects in the point particle EFT must also 
involve a discussion of absorptive effects, which for spin- 
less black hole binaries in the extreme mass limit are 
known to arise first at order Such effects are 

also encoded in point particle operators, although there 
are subtleties involved in the construction which we will 
discuss in a separate paper. 



VI. CONCLUSIONS 

In this paper, we have reformulated the NR expansion 
for binary star systems, an important class of signals for 
the gravitational wave detection experiments, in the lan- 
guage of effective field theory. The EFT approach is use- 
ful for the binary problem because it allows the different 
length scales that arise in the binary star to be treated 
independently. In particular, it is possible to separate 
the model dependent aspects of the problem which in- 
volve the internal structure of the binary constituents 
from the gravitational dynamics that governs the orbital 
evolution and radiation emission. Using the EFT for- 
mulation, it is easy to show that the divergences which 
arise at in the PN expansion can be attributed to the 
existence of new operators in the effective point particle 
description. However, these operators can be removed 
via a point transformation of the metric tensor and thus 
never contribute to physical quantities. This leads to the 
conclusion that there are no finite size effects at order v^. 
Practically, this means that whenever one encounters a 
log divergent integral at order in the potential, one 
may simply set it to zero. Its value cannot affect physical 
predictions. 

Beyond order v^, a new set of operators, which are 
quadratic in the curvature and whose role is to describe 
tidal deformations (induced moments), contribute at v^*^. 
The magnitude of such tidal forces is consistent with con- 
ventional Newtonian estimates of their size. However, in 
the EFT these conclusions follow from general princi- 
ples (dimensional analysis and general covariance) with- 
out the need to carry out any calculations. As will be 
discussed in a future publication, the calculations neces- 
sary to determine the coefficient of these operators (which 
encode the finite size effects) are relatively straightfor- 



ward, since they can be done for a compact object in 
isolation without the need to include superfluous infor- 
mation about the complete two-body problem. This is 
the beauty of working within an EFT formalism. 

Several aspects of the formalism remain to be worked 
out. Here we have not discussed logarithmic divergences 
that may arise in NRGR due to graphs with internal ra- 
diation graviton lines. Such divergences can be treated 
by the same renormalization and matching procedure de- 
scribed in this paper to deal with divergences in the cal- 
culation of potential terms. In fact, we believe that the 
explicit decoupling of the orbital scale r and the radiation 
scale r/v will make the evaluation of effects such as radia- 
tive tail terms, which is difficult within conventional PN 
approaches j^i^], simpler within our approach. Also, in 
our analysis so far, we have neglected the role of spins or 
higher rank tensors carried by the point particles. Such 
permanent moments encode the intrinsic, rather than the 
tidally induced, shape of the binary star constituents. It 
is important to extend the formalism, in particular the 
velocity power counting rules, to systematically include 
the effects of these additional worldline degrees of free- 
dom on the pattern of gravitational waves seen by the 
detector. We expect also that the inclusion of these dy- 
namical quantities will allow the construction of a vastly 
larger set of coordinate invariant worldline non-minimal 
couplings than those considered in this paper, and lead to 
an even richer structure of coupling constant RG flows. 
Finally it should be interesting to apply the formalism 
to work out the phenomenological signatures of theories 
of gravity that predict deviations from general relativ- 
ity, including simple extensions with additional fields or 
more speculative models that deviate from conventional 
gravity in the far infrared [ssl l36l| . 
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Note added: While this paper was being reposted, 
Ref. appeared which discusses the same field redefi- 
nitions for removing the operators / drR, J drv^^v'^ Rf^u 
mentioned above. 



APPENDIX A: POWER COUNTING L 

In this appendix we show that any diagram contribut- 
ing to Sf.ff[xa\ scales at most as L^, up to powers of 
the velocity. Consider a vacuum diagram with Ng gravi- 
ton self-interaction vertices from the Einstein-Hilbert ac- 
tion, insertions of the point particle worldline, and P 
graviton propagators. Writing the graviton kinetic term 
as ~ TOpj / d^x{dh)'^, each graviton self-interaction scales 
as m?pi , a coupling to the worldlines scales as the particle 
mass m, and each graviton propagator scales as 1/mp;. 
Up to powers of v, the only other scale that can appear 
in a Feynman diagram is the orbital distance r. Thus a 
typical diagram in the expansion of S'e/ / [xa] has a mag- 
nitude 



2Ng -2P N„ 



+2N„-2P 



(Al) 



since Seff[xa] is dimensionless. Using the virial theorem, 
nipf ~ v'^r/m, so up to powers of v, this equation im- 
plies that each diagram scales as {rnr)^^^^i'~^ ~ L^~^, 
where V = N^+Ng is the total number of vertices in each 
connected vacuum diagram that contributes to Seff[xa]- 
Starting with the trivial diagram with V — I, P — 
(corresponding to a free worldline) one can attach to it 
enough vertices and graviton lines to build up a diagram 
with arbitrary V, P. However, since the diagram obtained 
in this way must be connected, for each new vertex we 
must include at least one graviton line, so that V ~ P 
for this new graph is bounded by unity. From this we 
conclude that any diagram contributing to Sef / [xa] scales 
as L^~^, with integer £ >0. 



APPENDIX B: FEYNMAN RULES 

In this appendix we review the derivation of the Feyn- 
man rules used to obtain the results of this paper. As 
explained in the text, in our formalism all observables 
are derivable from either the quantity Seff[xa] defined 
by Eq. JBJ or from the effective Lagrangian for the radi- 
ation graviton mode, SNRGR[h,Xa]- Note that because 
there is no propagator for the "fields" a;^(T) (as far as 
the gravitons are concerned, these degrees of freedom are 
treated as a fixed background) we have the rule 

• iSeff[xa\ is the sum of Feynman diagrams that re- 
main connected when all particle worldlines are re- 
moved. 

Thus diagrams like the Newton exchange term in Fig.^t- 
are kept in the expansion, but diagrams such as a "box" 



diagram with a pair of worldlines interchanging two 
gravitons do not arise in the expansion of Seff[xa]- In 
fact, in our formalism the box diagram corresponds to 
the product of two single exchange terms, and is already 
accounted for by keeping just Fig. 

The procedure for integrating out modes with wave- 
lengths shorter than the radiation scale r/v involves com- 
puting diagrams with one or more external radiation 
gravitons and any number of internal potential gravitons. 
This gives rise to the quantity SNRGR[h, Xa] defined in 
Eq. 1)14(1 . The rule for calculating this is 

• A term with n powers of the radiation field ft.^^ 
in iSNRGR\h,Xa] is equal to the sum of Feynman 
diagrams with n external radiation graviton lines. 

Due to the multipole expansion, this quantity is guar- 
anteed to be of the general form SNRGR[h,'Ka] = 
J dti[/i(Xcm, t), Xa], as we explicitly saw in the deriva- 
tion of the low order terms in Seff[h,Xa]- 

We now give the Feynman rules for computing the di- 
agrams in the expansion of iSeff[xa\ and iSNRGR[h, Xa]. 



1. Propagators and external gravitons 

If we expand the metric perturbation h^^, = g^^ — rj^^ 
in terms of potential and radiation graviton modes as 



hf,^{x) = hf^i,(x) + 



(Bl) 



then the propagators for H\^^^ and h^j^i, can be read off 
the terms in the gravitational action that are quadratic 
in the fields. In order to get a well defined propagator 
one needs to fix the gauge, which is achieved by adding a 
term to the action that explicitly breaks diffeomorphisms 
This procedure is standard for gauge theories and can be 
found for instance in chapter 16 of Ref. ^J. For gravity 
this is explained in the reviews [2ll |2^ . Given the choice 
of gauge fixing term for potential gravitons used in the 
text, Eq. (|15l) the terms quadratic in the potential mode 
can be written as 

Sm = -\ j dx°^i/_k^,(x")k2T^''^"^7?k.^(x"), 

(B2) 

withT'^'''"^ = \■q^'°'■q'''^ + \r]^'f'r]''°'-\r]^"''q°''^ . The prop- 
agator is then defined as a formal inverse inverse to this 
equation 

,}^T^^'^;p^H^^Ax^)H^^^{x'^)) = 

(2^)353(k + q)^(:E0-x'0)/^;, (B3) 

with /'^"•"^ = ^T^'^^Tj-'f^ + ^Tj'^f^r]"" the identity on sym- 
metric two-index tensors. Inverting we find Eq. (|17|) 

(i/k^,(a;0)i/q„^(0)) = -(27r)^<53(k + q)j^5(x")P^,;„0. 

(B4) 



-r 



2. Vertices 



FIG. 11; Vertex associated with the Feynman rule in Eq. IIB9I 



The propagator for the radiation graviton also depends 
on the choice of gauge fixing term in the action. For 
instance, choosing a gauge fixing term of the form 



Sgf [h] = m%i J d'^xr^^'f^t, , (B5) 

with f ^ = <9q/i^ — ■^df^h, the propagator is given by 

{h^,i,{x)hat}{v)) = Df{x - y)Pfj.u;ai3, (B6) 
where the Feynman propagator D p (x) is given by 



Df{x) 



{2-kY + w 



(B7) 



In this equation, e — > 0+ defines a choice of contour for 
the integral. Given these results, the general rule for 
a given Feynman diagram is then: 

• An internal potential graviton line with label mo- 
mentum k and connecting points at times and 
x'^ corresponds to an insertion of Eq. ljB4|l . 



• An internal radiation line connecting two points x 
and y denotes a factor of Eq. ljB6|) . 



• For each external radiation line starting at a vertex 
located at a spacetime point x, one includes a factor 
of the radiation field hi_iv{x) and an integral J d'^x 
over X. 



They Feynman rules for general relativity contain an 
infinite number of interaction vertices. This is because 
the Einstein- Hilbert Lagrangian ^/gR is an infinite series 
in the graviton h^,y, so that the set of terms in this expan- 
sion with n factors of /i^;^ gives rise to an n-graviton ver- 
tex in Feynman diagrams. In our theory, the expansion 
of the point particle couplings also generates an infinite 
set of vertices that can be used to build up Feynman di- 
agrams. The method for deriving the vertex rules associ- 
ated with these interactions is straightforward, although 
in the case of graviton self-interactions arising from the 
gravitational action the procedure gets messy due to the 
large number of terms at a given order in /i^j^. In this 
section we will illustrate the procedure by deriving some 
of the vertex rules used in the text. The derivation of 
more complicated Feynman rules is directly analogous. 

Consider first Feynman rules arising from insertions of 
point particle interactions. The simplest non-trivial term 
is the vertex that arises in the calculation of the two-body 
Newton interaction. Expanding the — m J dr term in the 
action, there is a term of the form 



5. 



pp 



2mpi 



dx" [ e^''-(-")i/koo(x") + ---. (B8) 



This leads to a Feynman rule for the vertex shown in 
Fig. 11 



dx^y^e^'^-f-^'^.orM. (B9) 



Then for instance the diagram of Fig. [5^ is given by 



J 



Fig.|3l: 



imi 
2mpi 



dx° I e^i'-^if^ ^rj,,oV.o 



im2 
2mpi 



dx" 



{K''{x^)Hf{x'^)), (BIO) 



r 



which reduces to Eq. 1)22(1 . Likewise, the term in the 
action 

S,, = J dx" J e'''-(-")v,(x")iIko.(a;°) + • • • , 

(Bll) 



leads to the Feynman rule 



which can be used to calculate Fig.^Jj. The point particle 
action has also terms that are non-linear in the graviton 



and lead to vertices with more graviton lines. For in- 
stance, the vertex in Fig.jSja involving two graviton lines 
arises from a term in the action 



8Tn% 



dx° / e^('i+'')-(-'')i/koo(x°)i/qOo(x°). 



(B13) 

The Feynman rule for this term contains two identical 
terms, which arise from the two independent ways of con- 
tracting the two fields in the interaction with the external 
graviton lines (in other words, of pairing the two graviton 
lines with index pairs ^v, a[3 with the two fields in the 



22 



interaction) 



dx° / e'^^+'^y^^-°^f^o^.Vo. x 7joc.mf^- (B14) 



k,q 



Using this vertex, we can then calculate the graph in 
Fig. Including the proper symmetry factors for the 
diagram (see for instance chapter 4 of Ref. [33) leads to 

Eq. jsni)- 

Besides the interactions vertices involving point par- 
ticle sources, we also need to calculate the n-graviton 
vertices. Consider for instance the 3-point vertex for ra- 
diation gravitons. Due to the large number of terms each 
with a complicated tensor index structure, we do not dis- 
play the full Feynman rule here (it is more efficient to do 



the calculations using an algebra manipulation computer 
package). Consider, however one of the terms that gives 
rise to the 3-graviton Feynman rule. In the gauge of 
Eq. (|B5|) all terms in this vertex derive from the cubic 
terms in the expansion of y/gR. For instance, using 



2m% 



(B15) 



there is a term in the Einstein-Hilbert Lagrangian of the 
form 



2ml, I d'^x^R = — f d'^x - -h^ah'^'^ R'-^^ + — f 

J 2mpi J [ 2 \ mpiJ 



(B16) 



r 



where R^,}} = ^d^^dah" + ^d^^dah" — \d^h^^ — ^d^^d^h, term leads to a contribution to the 3-graviton vertex of 
and = -q^^R^ilu- In coordinate space this particular ^^'^ ^^^^ 

I 



mpi 



-T " / 



2 Ml''! a^A'2i'2 f3 

I 



perms. 



(B17) 



The tensor Ip,^af3 has been defined below Eq. HB3|) . In all permutations of the terms shown in the above expres- 

this vertex, each incoming graviton is labeled 1, 2, or sion under the interchange of the graviton labels. It is 

3. The partial derivatives in this particular term act on also useful to work in momentum space, in which case 

the propagator of the graviton labeled 3 coming into the the expression for this piece of the 3-graviton Feynman 



vertex. To get the full Feynman rule from the term in 
Eq. IjBlGjl . Bose statistics also demands that we include 



rule becomes (taking all graviton momenta to be flowing 
into the vertex) 



mpi 



{2^Y5''{ki +k2 + h) 



-I "7 



1^ -I-Ia aT'^P \ T ( hP h"^ - nP" h'^\ 



^ -^mvi „-^P2'^2aP I „^ M3f3'>'3 1^3 a ^ r, ^ A'31'3'^3 l^3a r, ^ fJ-s'^s'^S r, P- 



a '^3 H 



+ perms. (B18) 



Feynman rules for potential self-couplings, or for potential-radiation couplings are derived in the same way. 
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